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Abstract 

The (A^3^) decay has been measured with stopped positive kaons 

at the KEK 12 GeV proton synchrotron. A sample containing 125 events was 
obtained. The partial branching ratio Br{Kfj, 3 j, E.y > 30 MeV, 7+7 > 20°) was 
found to be [2.4 ± 0.5{stat) ± 0.6(syst)] x 10“®, which is in good agreement with 
theoretical predictions. 
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Semi-leptonic radiative decays of K-mesons, K nlwy offer a good 

testing ground of hadron structure models making use of low-energy effective 
Lagrangians inspired by chiral perturbation theory (ChPT). It is expected 
that branching ratio measurements and decay spectra with a single pion in the 
hnal state provide simple but good constraints on the models. The radiative 
decays of mesons usually consist of an internal bremsstrahlung (IB) process 
and a hadron-structure-dependent direct emission (DE) process. While the 
IB process is dominant in the decays with electrons in the hnal state such 
as K ^ TTCz/y decays, one expects a signihcant DE contribution when 

there is a muon in the hnal state, K Tr/iz/y because of the larger 

lepton mass. The relative size of the DE ehects can be calculated in strong 
interaction models. Following the early estimates [1,2,3], based on current 
algebra, calculations in the framework of the ChPT theory have been done [4]. 

branching ratios of neutral kaons with electrons and muons in the hnal 
state have been reported in the literature with branching fractions of 
-^°37 of 3.5 X 10“^ and 5.5 x 10“"^, respectively [5]. For the charged kaons, a 

decay branching ratio of 2.65 x 10“^ has been measured [5], and results 
of the hrst measurement of the decay using an in-hight K~ beam has 
recently been reported [6]. In this paper, we present a new measurement of the 
7r°/i+z/7 decay using a stopped beam along with detailed 

Monte Carlo simulations, which enabled us to determine the branching 
ratio. 

The experiment was performed at the KEK 12-GeV proton synchrotron. The 
detector was basically the E246 setup [7], which had the 12-sector toroidal 
spectrometer and the ancillary detector assemblies such as the photon calorime¬ 
ter and the particle tracking system. Since the system was built primarily 
for the purpose of a high precision test of time reversal invariance in the 
7r°/i+z/ decay [7], an elaborate simulation program based on 

GEANT3 [8] has been developed. Details of the setup are well documented in 
Ref. [9]. In addition to the T-violation search, spectroscopic studies for various 
decay channels have also been successfully performed using the same detector 
system [10,11,12]. 

A separated 660-MeV/c beam was stopped in an active target system. The 

Kfj.S'y events were identihed by analyzing the /i"*" momentum with the spectrom¬ 
eter and detecting three photons in the CsI(Tl) calorimeter. The momentum 
vectors of the charged particles were determined by reconstructing their trajec¬ 
tories in the spectrometer using multi-wire proportional chambers (MWPCs). 
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The fi'^s were selected by determining the squared mass from a time- 

of-flight measurement. The photon energy and hit position were obtained, 
respectively, by summing the energy deposits and taking the energy-weighted 
centroid of the CsI(Tl) crystals sharing a shower. The analysis procedures of 
the present work for the charged particle tracking, TOF measurement, and 
photon energy and hit position determinations are the same as those of the 
previous K.„ 2 'y study (see Ref. [11] for details). 

Specihc cuts for the selection are described below. The charged particle 
momentum corrected for the energy loss in the target (T)j+) was imposed 
to be P^+ < 170 MeV/c. Events from 7r~^ decays in-flight and scattering of 
the charged particle from the magnet pole faces were eliminated by requiring 
the particle track to be consistent with the hit position in the ring counters 
surrounding the active target system [9]. The selection criterion for muons was 
8000 < MpQp < 14500 MeV^/c^, as shown in Fig. 1. Events with three photon 
clusters in the calorimeter were selected: two as coming from 7 r° ^ 7172 and 
one being a radiative photon ( 73 ). Since there are three possible combinations 
to form a 7 r° from three photons, a quantity was introduced to hnd the 
correct pairing, 

Q2 = _ Mf/al^ + (cos0jl^4 - COS0j^4 “ (1) 

where M^-o is the invariant mass of the selected pair and 9^+^2, is the opening 
angle between the and 73 . The superscripts MEA and CAL stand for the 
measured angle and the angle calculated from other observables by assuming 
the kinematics. The pair with the minimum {=Qmin) among the 

three possible combinations was adopted as the correct pairing. The a {<Jm, 
(Tq.) and offset values (M, a) in each terms are aM = 10.92 MeV/c^, cTq, = 
0.273, M = 118.3 MeV/c^, and a = 0.265. The choice of the parameters were 
determined to obtain the highest probability for the correct pairing by using 
the simulation data. The correct pairing probability was estimated to be 69% 
from the Monte Carlo simulation. Further, since most of background events do 
not satisfy the iF^ 3 .y kinematics, the cut of Q%n <1-5 reduced the background 
contaminations. An additional cut condition, cos[ 6 *.^.^]min < 0.45, was applied 
to reject events with a photon split into multiple clusters, where [^^^jmin is 
the minimum opening angle of photons among the three combinations. The 
above conditions were sufficient to select events. From this analysis, a 

sample with 565 events was extracted. The spectra are shown in Fig. 2. The 
black(sohd) histograms are the data comprising of the i 7 ^ 3 .y events and the 
background events to be discussed below. 

There are three major background components, K~^ 7r~^7r^7r^ {K^s), 

TT+TfOy (^KT^ 2 'y), and 77 ^ 3 . The former two could imitate if the pion decays 
in flight and the three photons hit the calorimeter. Also, 77^3 with an acciden¬ 
tal photon could contribute to 77 ^ 3 .y. The 77,^3 and 77,^27 contaminations were 
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Fig. 1. Correlation plot of and The i^Tra and K-,^ 2 'y events which were used 

to calculate the branching ratio and the background fractions are also seen. 

estimated using a Monte Carlo simulation. The simulation data were analyzed 
in the same manner as the experimental data, yielding surviving background 
fractions. In order to determine these fractions, the results of careful evalua¬ 
tions, carried out in the previous i^,r 27 study [11], were used. These and 
i^ 7 r 27 events can be seen in the scatter plot in Fig. 1. The numbers 

of the iCjrs and ^^,^27 events were calculated from those of the experimental 
and iF 7 r 27 events by using acceptance ratios as. 


y{KT) 






( 2 a) 

( 2 b) 


where Y{X) is the yield of decay channel X and f2(X) is the detector ac¬ 
ceptance determined by the simulation. BG and NM stand for the selection 
conditions of the present background evaluation and the previous normal i^ 7 r 27 
study [11], respectively. Potential systematic errors from the uncertainty of the 
calculated acceptances can be reduced by taking the ratio of the acceptances. 
Substituting the geometrical acceptance of our setup, = 9.91 x 10“’^, 

= 9.45 X 10-^ x lO”®, and = 2.40 x lO'^ 

and the measured yields Y (iF^^)=32105 and Y{K^^) = 4434, Y{K^^) and 
y (-^ 7 ^ 27 ) were determined to be 337 and 67, respectively, corresponding to 60% 
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Fig. 2. Black(solid) histograms are experimental data extracted with the present 
selection conditions for (a) (b) Q^, (c)M|.Qp, (d) E-y^, (e) cosO^+y^, and (f) 

Eg. Red (dot-dashed), blue(dashed), and green(dot) histograms are backgrounds due 
to RjtS, ET^ 2 'y, and decays, respectively. 


and 12% of the total events. The contribntion from was stndied by ac¬ 
cepting accidental events in the CsI(Tl) TDC data and its nnmber was fonnd 
to be 36. The backgronnd spectra are shown in Fig 2 as the red(dot-dashed), 
blne(dashed), and green(dot) histograms for KT, 2 'y, and respectively. 
Fig. 3 shows the backgronnd-snbtracted distribntions compared with the sim- 
nlation described in the following. The nnmber of K^^y events was fonnd to 
be 125 ± 25. 
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Fig. 3. spectra obtained by subtracting backgrounds (dots): (a) P^+, (b) , 

(c)M|.qp, (d) £^^ 3 , (e) cos0^+.y3, and (f) Eg. The histograms are the Monte Carlo 
simulations. 


In the simulation, the data were generated according to the matrix 

elements given in Ref. [1]. Here, only the IB process was taken into ac¬ 
count because the present result is not sensitive enough to study the DE 
contribution (see below). The detector acceptance for in fhe region of 
E.yg > 30 MeV, 6 '^+.y 3 > 20° was determined to be r 2 (iC^ 3 ..y) = 1.52 x 10“^. The 
reduced values of hts between the experimental data and the simulation 
are (a)1.33, (b)1.09, (c)0.80, (d)1.24, (e)0.94, and (f)1.05 in each spectrum, 
lending credence to our claim that these are indeed events. 
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The partial branching ratio of in the region of > 30 MeV, 6*^+^3 > 20° 
can be derived from separately comparing with the i^,r 27 and branching 
ratio by means of the yield and acceptance ratios as, 






Y{K^^^) n{K,s,) 
Y{K,s,) 


Bt i^K■j^2'y)) 


BriK^z). 


(3a) 

(3b) 


where Br{X) is the branching ratio for the process X. Making use of i?r(if^ 27 ) = 
2.61 X10“"^ predicted by theoretical calculation [13] and i?r(i^^ 3 ) = 1.73 x 10“^ 
by the Particle Data Group [5], the branching ratio was determined to 
be Br{Kf^ 3 .y) = (2.4±0.5) x 10“^ from iP,r 27 and Br{K^ 3 ^) = (2.8±0.6) x 10“^ 
from 77 ^ 3 , which are consistent with each other. Since only tt+s in the end¬ 
point region of the 77^3 decays were accepted by the toroidal spectrometer 
and the spectrometer acceptance was strongly affected by the 77^3 form fac¬ 
tors assumed in the simulation, we adopted 7 ?r( 77 ^ 3 .y) from Eq.(3a) and used 
Eq.(3b) only for a consistency check of the deduced results. 


The major systematic errors in the determination of the 77^37 branching ratio 
come from uncertainty of the 77^3 and 77^27 background fractions. Since we 
rely on the Monte Carlo simulation for the background estimation, imperfect 
reproducibility of the experimental conditions introduces systematic errors. 
We carefully evaluated the effects from the TOE measurement and the charged 
particle tracking, which would distort the 77 ^ 3 .y, 77 ^ 3 , 77^27 spectra and, as a 
consequence, introduce the systematic errors. These errors were estimated by 
varying the inputs of the TOE resolution and the MWPC spatial resolution in 
the simulation within their maximum likely uncertainties (see below). In order 
to estimate the other systematic effects, the 77^3 and 77^27 fractions could be 
controlled by requiring an extra cut at the cost of good 77 ^ 3 .y events as follows. 

Since events with tt’*' decays in the spectrometer (DIS) are seen as a tail 
under the tt’*' peak, as shown in Fig. 2(c), a wrong estimate of the TOE 
resolution could introduce a systematic error. This effect was estimated by 
changing the TOE resolution assumed in the simulation from utof = 270 
ps, which was experimentally determined using the Kes and 77^3 decays, to 
o'TOF = 270 ±20 ps. Its contribution to the branching ratio error was found to 
be ABr{K^ 3 ^) = 0.1 x 10 “®. Also, the muon mass cut point dependence was 
investigated to study its contribution to the DIS component and found to be 
very small. We note that an asymmetric structure due to the DIS contribution 
shown in Fig. 2(c) was removed by subtracting the 77,^3 and 77^27 backgrounds. 

Since most of DIS events were rejected by the consistency cut between the 
hit position in the ring counter and the charged particle track, the MWPC 
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spatial resolution for the input to the simulation affects the estimation of the 
background fraction. In the simulation, we varied the resolutions by ±10% 
from the experimental values to estimate this effect. The background fraction 
was obtained by repeating the same analysis. We arrived at ABr^K^j^^^) = 
0.2 X 10“^ as the contribution to systematic error from particle tracking. 

The background contribution is also influenced by the total energy sum of 
the /i’*' and the photons dehned as Eg = E^+ ± E^i because the Eg 
distribution has a specihc structure corresponding to the decay channel, as 
shown in Fig. 2(f). Selecting events with the conditions of Eg > 400 MeV 
and Eg < 400 MeV, the branching ratio was determined to be (2.9 ± 
0.6) X 10“® and (2.0 ± 0.7) x 10“^, respectively. The background fractions 
using these cut conditions are summarized in Table 1. It is noted that the 
Eg variation affects the background contributions signihcantly but does not 
change the branching ratio which depends on the ratio of events and 

acceptances. This feature indicates the correct estimation of the and K-,^ 2 'y 
background fractions and the detector acceptances. Although these parameters 
are consistent within errors, the parameter shifts, ABr{K^ 2 ,'y) = 0.5 x 10“^, 
were treated as a systematic error due to an acceptance uncertainty for the 
background events. 

The 7 mispairing effect was studied by varying the pairing probability. In addi¬ 
tion to the standard selection conditions, an extra cut, cos 6 *y,i.y 2 > —0.5, 
was imposed, where 0^^^^ is opening angle between 71 and 72 , and the pairing 
probability became 81% with a loss of good events,. The Br{K^ 3 ^) was 
found to be consistent within 0.2 x 10 “^ which was regarded as the system¬ 
atic error. The statistical fluctuation of the backgrounds from with an 
accidental photon was treated as a systematic error of 0.05 x 10“^. The total 
systematic error was evaluated by adding all these contributions in quadra¬ 
ture, as shown in Table 2. The most dominant effect is due to the error of the 
background fractions from the uncertainty of the detector acceptance. 

In conclusion, we have performed a measurement of the radiative de¬ 

cay using stopped positive kaons. The data sample of 125 events was ob¬ 
tained by subtracting the 77 ,^ 3 , -^^ 27 , and backgrounds. The hnal result 
for the branching ratio is Br{K^ 3 ^, E^^ > 30MeV, 6 *yj+.y 3 > 20°) = [2.4 ± 
0.5{stat) ± 0 . 6 (s|/st)] x 10 “®, by normalizing to the theoretical value of inter¬ 
nal bremsstrahlung in the iF,r 27 decay. This result is in agreement with the 
prediction from chiral perturbation theory [4]: E^^ > 30MeV, > 

20°) = 2.0 X 10“®. The current experimental result is not yet sensitive enough 
to study the contribution from hadron structure effects, which is predicted 
to be at the level of a few % of the rate [4], and/or to search for T-violating 
triple correlations such as Pt = p-y- xp/o) [1,14]. It is conceivable that such 
experiments may be attempted at the new J-PARC facility currently under 
construction [15]. 



This work has been supported in Japan by a Grant-in-Aid from the Ministry of 
Education, Culture, Sports, Science and Technology, and by JSPS; in Russia 
by the Ministry of Science and Technology, and by the Russian Foundation for 
Basic Research; in Canada by NSERC and IPP, and by the TRIUMF infras¬ 
tructure support provided under its NRC contribution. The authors gratefully 
acknowledge the excellent support received from the KEK staff. 


9 



References 


[1] E. Fischbach and J. Smith, Phys. Rev. 184 1645 (1969). 

[2] H. W. Fearing, E Fischbach, and J. Smith, Phys. Rev. D 2 542 (1970). 

[3] H. W. Fearing, E Fischbach, and J. Smith, Phys. Rev. Lett. 24 189 (1970). 

[4] J. Bijinens, G. Ecker, and J. Gasser, Nucl. Phys. B 396 81 (1993). 

[5] Particle Data Group, Phys. Lett. B592 1 (2004). 

[6] O. Tchikilev et al, hep-ex/0506023 (2005). 

[7] M. Abe et al., Phys. Rev. Lett. 83 4253 (1999); M. Abe et al. Phys. Rev. Lett. 
93 131601 (2004). 

[8] R. Brun et al., CERN Program Library Long Writeup W5013, CERN 
Applications Software Group (1993). 

[9] J.A. Macdonald et al., Nucl. Instr. and Meth. A506 60 (2003). 

[10] S. Shimizu et al., Phys. Lett. B495 33 (2000); A.S. Levchenko et al., Phys. At. 
Nucl. 65, 2232 (2002); K. Horie et al., Phys. Lett. B513 311 (2001); Y.-H. Shin 
et al., Eur. Phys. J. C12 627 (2000). 

[11] M.A. Aliev et al, Phys. Lett. B554 7 (2003). 

[12] S. Shimizu et al., Phys. Rev. D70 037101 (2004). 

[13] G. D’Ambrosio, M. Miragliuolo, and P. Santorelli, in Da^ne Physics Handbook, 
edited by L. Maiani, G. Pancheri, and N. Paver (Laboratori Nazionali di 
Frascati, Frascati, 1992). 

[14] V. V. Braguta, A. A. Likhoded, and A. E. Chalov, Phys. Rev. D 68 094008 
(2003). 

[15] “The Joint Project for High-Intensity Proton Accelerators”, KEK Report 99-4 
(1999); http://j-parc.jp/ 


10 



List of Figures 


1 

Cnrrdalion niot of M^r^r, aud P... The K^-:, a,nd events whieh were used to ealeniate the 


2 

Ria.ektsoiid'l histograms are exnerimental data extracted with the oresent seleetion eondition 


3 

soectra obtained bv subtracting backgrounds fdotsk fa) P,,+ . ibi O^. fc)M'^r^p. fdi 


11 











Table 1 

Variation of the and K-,^ 2 '^ background fractions with the Eg cuts. The Y 
and Br{Kfj^^^) values are also shown. 


^.(MeV) 

K^3 

K'7r2'y 


Br{K^3^)xl0^ 

no cut 

60% 

12 % 

125 

2.4 ±0.5 

Eg > 400 

28% 

27% 

68 

2.9 ±0.6 

Eg < 400 

74% 

5% 

57 

2.0 ±0.7 


Table 2 

Summary of major systematic errors. All items are added in quadrature to get the 
total. _ 


Error source 

Uncertainty of ilr(Ar^ 3 .y) xlO^ 

TOE resolution 

0.1 

MWPC resolution 

0.2 

Detector acceptance 

0.5 

7 mispairing 

0.2 

Accidental background 

0.05 

Total systematic error 

0.6 
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